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We analyze the spectral properties of a phenomenological model for a weakly doped two-dimensional anti- 
ferromagnet, in which the carriers move within one of the two sublattices where they were introduced. Such 
a constraint results in the free carrier spectra with the maxima at k = (±7r/2, ±7r/2) observed in some 
cuprates. We consider the spectral properties of the model by taking into account fluctuations of the spins 
in the antiferromagnetic background. We show that such fluctuations lead to a non-pole-like structure of the 
single-hole Green's function and these fluctuations can be responsible for some anomalous "strange metal" 
properties of underdoped cuprates in the nonsuperconducting regime. 
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1. Introduction 

After more than twenty years from its discovery, the problem of high-temperature supercon- 
ductivity (HTSC) remains unsolved. Nevertheless, there are some facts about HTSCs, which are 
generally accepted by the scientific community. In particular, it is well known that these mate- 
rials transform from antiferromagnetic insulators into superconductors with carrier doping. The 
superconductivity in most of the cuprates mainly takes place in the CuC"2 layers, and the other 
inter-layer atoms supply the carriers for these layers and play a role of the carrier scatterers. It 
is believed that the presence of the antiferromagnetic background strongly affects the physical 
behavior of the weakly doped materials. In particular, this behavior can be defined by strong hole 
or electron correlations [HI [2]. Moreover, the correlations can be responsible for the <i-wave super- 
conductivity in many cuprates [OS]. In fact, since the isotope effect in optimally doped HTSCs 
is rather weak, it suggests that the electron-phonon coupling is not the main source of supercon- 
ductivity in these materials, though the role of phonons and, in particular, the interplay of strong 
correlations and phonon coupling in HTSCs is currently an active area of research (see, for exam- 
ple, Ref. [[5] and references therein). Thus, it is believed by a significant part of the researchers 
that the phenomenon of HTSC can be explained by using a strongly correlated model, in which the 
superconducting pairing with the anisotropic order parameter is caused by an antiferromagnetic 
spin-wave coupling. Probably, the most popular models for this scenario are the two-dimensional 
Hubbard model and its approximation in the case of strong correlations, the iJ-model [HIE] (see, 
e.g., Refs. [HE [HIED [T2]). 

Unfortunately, these models cannot be solved exactly in the two-dimensional case, so it is dif- 
ficult to make a firm conclusion whether they can be considered as realistic models of HTSCs. As 
an alternative approach, one can consider a simplified phenomenological model of cuprates with 
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strong correlations, which takes into account their main properties and can be solved exactly. Typ- 
ical unusual and important properties of cuprates include a universal dome shape of the critical 
temperature - doping curve, a different from the BCS theory ratio 2A(T = 0)/T™ ax ~ 5.5 at 
optimal doping, different critical densities for the superconducting gap in the nodal and antin- 
odal directions. Probably, the most unusual phenomenon in HTSC is the pseudogap phase in 
the underdoped regime. In this phase, the materials demonstrate very unusual properties of a 
"strange metal", like an anomalous temperature dependence of the resistivity etc, which are dif- 
ferent from a Fermi-liquid behavior. There is no general agreement in the HTSCs community on 
the origin of this phenomena, however a recent improvement of the angle-resolved photoemission 
spectroscopy (ARPES) technique gives a hope that the mystery of the "strange metal" phase in 
underdoped cuprates will be resolved soon. Another unusual feature of some cuprates is their free 
carrier spectrum. For example, the free spectrum of Sr2Cu02Cl2 and some other materials in the 
insulating phase has the maxima at the momenta k = (±7r/2, ±7r/2) [[T3]. This fact suggests that 
the main hopping processes in these systems take place between the next nearest neighbor (NNN) 
and next NNN sites. In other words, they correspond to an inter-sublattice carrier motion. Since 
the NNN and next NNN hopping parameters are too small for the oxygen sublattice of the CuC>2 
planes, it is difficult to believe that the holes move within the oxygen sites, as it is assumed by many 
researchers [ [T3] . It is also difficult to assume that the free carrier spectrum corresponds to the 
copper site NNN and next NNN hoping of the Zhang-Rice singlet, formed by an oxygen hole and 
by one of the copper ions. In fact, the Zhang-Rice singlet states, which move within their magnetic 
sublattice, are unstable due to the hole frustration with respect to the choice of the axis of the spin 
quantization and some other reasons [ 15, 16J. In order to describe some of the physical properties 
of the cuprates with the maximum of the free carrier spectrum at k = (±tt/2, ±7r/2), there was 
proposed a model, in which the holes occupy and move within the copper ion antiferromagnetic 
sublattice where they were born. [ 115) It was assumed that the hopping takes place between the 
NNN and next NNN sites, and the superconducting hole-hole attraction is due to the minimization 
of the energy of the system, when two holes occupy nearest sites. In this case, the minimal number 
of the antiferromagnetic bonds is broken [ [TTl, [TSl, [T9l [20l [2TJ and similarly to the Hubbard and the 
tJ-model cases, the pairing takes places predominantly in the cZ-wave channel [121]. 

We have already studied some of the properties of cuprates in the superconducting and pseudo- 
gap phases in the framework of this model. Namely, in Ref. [[IB] , we have analyzed the temperature- 
hole phase diagram of the model in the case of low doping. In this paper, we have obtained the 
doping dependence of the superconducting critical temperature T c by solving a system of coupled 
equations for the Green's function for the Hubbard operators within a generalized mean-field ap- 
proximation [ 22 . It was shown that superconductivity in the model arises at finite doping and 
T c grows with doping in the underdoped regime. We have also shown that there is an additional 
pseudogap phase at temperatures above T c and below another critical temperature To, which also 
grows with doping increasing in the case of low carrier density (for a schematic picture, see Fig. [I]). 
Namely, according to the Emery-Kivelson scenario [ [23] , the Cooper pairs start to form below 
the temperature To, which is associated with the superconducting mean-field critical temperature 
T^ F . In the two-dimensional case the pairs are disordered (the order parameter phases are ex- 
ponentially ordered) above the condensation temperature T c . This temperature can be associated 
with the Berezinskii-Kosterlitz-Thouless (BKT) temperature, below which the phases of the super- 
conducting order parameter are algebraically ordered. This is the only possible critical temperature 
in the two-dimensional case. The existence of the superconducting pseudogap phase with unusual 
properties at T c < T < Tq was confirmed in some cuprates, where a strong Nernst effect was 
observed (see, for example Refs. [[Ml US], and a theoretical paper [[55]). However, it is well-known 
that the anomalous properties of cuprates in the underdoped regime take place up to the tempera- 
tures much higher than Tq. In particular, an anomalous pseudogap in the one-hole density of states 
is observed below a temperature T* ^> T , which is called the pseudogap critical temperature. We 
present a schematic phase diagram of cuprates in Fig. [1] where we distinguish two regions in the 
pseudogap phase: the superconducting lower pseudogap (LPG) phase and the other, which we call 
the upper pseudogap (UPG) phase. The UPG critical temperature T* is a decreasing function of 



2 



Spin fluctuations and the "strange metal" behavior 




0.0 0.1 0.2 0.3 



carrier density 



Figure 1. A typical temperature-carrier density phase diagram of the hole-doped cuprates. It 
consists of the "normal", superconducting (SC), antiferromagnetic (AF), and two pseudogap 
phase regions: the lower pseudogap (LPG) phase, in which the behavior is affected by both SC 
and AF spin phase fluctuations, and the upper pseudogap (UPG) phase, in which only the AF 
spin phase fluctuations are present. 



doping. There are some experimental evidences that T* goes below T c in the overdoped regime and 
approaches zero at doping S ~ 0.19. It is believed by many researchers that the physical properties 
in the region T < T < T* are defined by nonsuperconducting processes (for over-view, see Ref. [ 
[21]). One of the popular explanations is based on the idea of the spin singlet formations in a doped 
two-dimensional antifcrromagnet, which corresponds to the resonant valence bond model. Unfor- 
tunately, there are no crucial experimental results which confirm the existence of such a state so 
far. Also, it is known that a non-Fermi liquid temperature dependence of the conductivity in the 
underdoped cuprates takes place up to temperatures of order 3000K, therefore, as it was suggested 
by Phil Anderson [ 27 , it shows that probably nonsuperconducting effects are responsible for this 
phase. Namely, he suggested that the unusual behavior of the spectral function can be explained by 
a renormalization of the quasiparticle Green's function due to a Gutzwiller projection in a strongly 
correlated model, which leads to an additional time-dependence of the Green's function and a 
non-pole-like (cut-like) form of this function in the frequency representation, which corresponds to 
a non-Fermi-liquid case. 

In paper [ [55] , we have studied some of the spectral properties of the model proposed in Ref. [ 
115] by taking into account fluctuations of the phases of the superconducting order parameter 
and the phases of the spins in the antiferromagnetic background. In particular, we have shown 
that the growth of the Fermi arcs with temperature in underdoped cuprates can be qualitatively 
explained within the model by taking into account fluctuations of the superconducting d-wave order 
parameter. In this paper, we analyze the effect of the superconducting and spin fluctuations on 
the structure of one-hole Green's function and the consequent anomalous behavior of the cuprates 
in the pseudogap phase. In particular, we show that a cut-like structure of the Green's function , 
qualitatively similar to the one obtained in Ref. [ 127] . can be obtained by taking into account these 
fluctuations. 

The paper is organized as follows. The model and the main equations in the mean- field case 
are presented in Section [2] In Section [3] we extend the problem by taking into account fluctuations 
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of the spins on the antiferromagnetic sublattices and estimate the doping dependence of the UPG 
critical temperature T*. The results for the spectral function and the density of states at different 
values of temperature are presented in Sections [4] and [5j correspondingly. In addition, in Section 
[5] we show how the anomalous frequency dependence of the conductivity can be obtained from 
the cut-like Green's function. A summary, a discussion of the results and conclusions are given in 
Section [6] 

2. Model 

As it was shown in Ref. [ [15] , the effective model for the holes in some of the weakly doped 
cuprates can be written as: 

H = (e d - M ) E ^ ~ \ E** m cos QAFM{ " " m) X?"XW J £ X?X% p . (1) 

11 11 . in n,p— a,b 

In the last equation, X 2,2 , X^ 2 ' 2 Xn' 1 ' 2 are the Hubbard operators for the hole number, the hole 
annihilation and the hole creation on the site n. The first three terms in the Hamiltonian Eq. (TTJ) 
describe the local hole energy and the NN, and NNN hopping processes, where Ed and \i are the 
hole on-site energy and the chemical potential. The magnetic structure vectors Qafm are equal 
to (±7r, ±7r), which corresponds to the antiferromagnetic case. We use the local spin coordinates 
for the Hubbard operators. The noninteracting part of the hole Hamiltonian describes holes, which 
move within their sublattices. The free hole dispersion relation is 

e(k) = Ed — 4i2 cos k x cos k y — 2*3 (cos 2k x + cos 2k y ) — fi. (2) 

The effective hole-hole attraction in the system is described by the last term in Eq. (fTJ). In fact, 
the doped holes introduced on the antiferromagnetic lattice, which move within the sublattice they 
were introduced, lead to a minimal increasing of the energy of the system when they sit on the 
nearest sites [ 117) , since in this case the minimal number of the antiferromagnetic couplings J 
between the nearest site spins is broken. In this case, two doped holes will always try to occupy 
NN sites, which results in the effective attraction described by the last term of the Hamiltonian 
Eq. (fTJ. In our calculations, we use the length units such that the lattice constant is equal to one, 
a = 1, and choose the energy parameter Ed to be equal to 4*2 + 4*3. In this case, the free hole 
energy is equal to zero at k = in the limit of low doping (/1 — > 0). In order to find dependence of 
the physical properties of the system on hole concentration 5, we shall use the following equation, 
which defines S in terms of the Hubbard particle number operator: 

S = J2(X 2 n > 2 ). (3) 

n 

Recently, we have studied the spectral properties of the model my taking into account super- 
conducting fluctuations [I28j. In this paper, we use a similar formalism to study the effect of the 
antiferromagnetic background spin fluctuations on the anomalous spectral and some other proper- 
ties of the system in the pseudogap phase at T > T c . Despite we are interested mainly in the spin 
fluctuation effects, we shall consider the general case assuming that the superconducting pairing 
can also take place, which corresponds to the temperature interval T c < T < Tq. At higher tem- 
peratures, i.e. at To < T < T*, we shall use the same equations by putting the superconducting 
gap to be equal to zero in our calculations. 

Then, in order to study the properties of the system described by the Hamiltonian Eq. (fTJ, it 
is convenient to introduce generalized Nambu-Hubbard hole operators 

*n(*) = f Sll® ) ■ = (x^(t),X^(t)) , (4) 
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where n is the lattice site and t is time, and to calculate the time-ordered Green's function 

G nm (t, f) - -i(T(* n (i)*L(*0))- (5) 
This function satisfies the following equation: 

ij t G nm (t, = 5(t - t')S nm i + (T[* n (t),H}*l(f)), (6) 

where / is a diagonal 2 x 2-matrix with the nonzero elements equal to (xi /21/2 (t) + X% 2 {t)). In 
the case of low doping, / ~ 1. In order to solve Eq. ([6]), it is convenient to approximate its last 
term by a generalized mean-field theory expression: 

(T[* n , fl]*L)(w) - E ^nlGlm(^), (7) 
1 

where 

E nm = ({[*„, (8) 

is the energy matrix (for details see, e.g., Ref. [ 122]). In this approximation, one neglects the 
dynamical corrections to the self-energy, which can be systematically taken into account. For 
example, in Ref. [[10] the authors considered these corrections in the case of the tJ-model by using 
a similar formalism. In order to find the Green's function, one needs to calculate the elements 
of the energy matrix E n i, which depend on different correlation functions, in particular on the 
superconducting gap function, which can be found by using the fluctuation-dissipation theorem. 
Assuming that the superconducting pairing takes place in the d-wave channel and introducing the 
relevant gap function: 

A d (k) = -4jE 7d (k) 7d (q)(x!' q /2 ^i/ 2) = Adld (k), (9) 
q 

where 7d(k) = cos(k x ) — cos(k y ) is a <i-wave structure factor, one can get the following approximate 
expression for the Green's function in the frequency-momentum representation: 

G ( w ' k ) = I 7U\ \ -a^ ' ( 10 ) 
lu + e(k)T z + iA(k)Ty 

where r y and f z are the Pauli matrices. 

In order to find the doping and temperature dependencies of the superconducting gap parameter 
Ad, one needs to derive and to solve the system of equations for and the chemical potential [i. 
These equations follow from the definitions Eqs. ([3]), ([9]) and the self-consistency conditions, which 
follow from the fluctuation-dissipation theorem: 



4J^(q)tanh vgM^M L__ 

(12) 



1 + tanh 



v^(qT+AH(q)\ e(q) 



2T I y?(^TAHM 



(see Ref. [[TB] for details). The solution of these equations at A^ = gives the doping depen- 
dence of the mean-field critical temperature T c Mi? , or the LPG critical temperature To. According 
to the Emery-Kivelson scenario, the real superconducting critical temperature T c < T c in the 
two-dimensional case corresponds to the BKT temperature, below which the phases of the order 
parameter become algebraically ordered (for over-review see, for example, Ref. [[21])- We studied 
the doping dependence of this temperature in Ref. [ 116) . In this paper, we shall mainly concentrate 
on the spectral properties of the model in the pseudogap phase at T > T c by taking into account 
superconducting fluctuations and spin fluctuations of the antiferromagnctic background. 
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3. Spin fluctuations 

In order to describe the UPG region, we assume that, in analogy with the superconducting 
order parameter fluctuations, the physics in the temperature range Tq < T < T* is governed 
by the fluctuations of the spin phases. Supposing that the antiferromagnetic copper oxide spin 
Hamiltonian is described by the XY-model, it is easy to obtain an analogous critical temperature 
for the spin subsystem. Therefore, we associate T* with the temperature of the BKT transition 
for spins. This temperature can be also estimated from the following equation: 

T* = ^J apin (»,T*), (13) 

where the spin stiffness J SV i n is the coefficient in front of the quadratic term of the spin phase 
gradients in the effective action for the spin phase ip n = Qafmh differences: 

n = *f*/ *r<p%<p + *±J d 2 r(W>) 2 . (14) 

It is equal to 

7 1 f ^ tanh ^ g " 2(k) + A ^ (k) 

J spin - ^ J (27r)2 tanh ^ 2T 

where e(k) = e(k) — 4J<5 + 4tiS(cosk x + cosk y ) (see [EH]). As it was mentioned above, doping 
leads to a gradual destruction of the antiferromagnetic order and to a hole hopping between NN 
sites, which belong to different sublattices. We take into account this process by adding the doping- 
dependent term 4i 1 5(cosfc c + cosfc y ) to the free spectrum Eq. ©. The doping dependence of T* 
can be found from Eq. (flU)) at = 0. In this case, Eq. (fT3"|) has the following form: 

4 J (2tt) 2 \2T* J y ' V ' 

The exact solution of Eq. (jTHJ) shows that T* is a decreasing function of doping. One can estimate 
the analytical dependence T*(S) from Eq. (j 1€3[) by taking into account that at low doping T* ^> 
|e(k)|. Therefore, in this case: 



T *^^lS (§M 2(k) - T * max{i - 1351 (17) 

where T* lax = / ^^(k)\ s=Q and /3 = 4^/2h J j0^e(k){cosk x + cosk v ). 



'■(k) 



vm~ 2 



(15) 



4. Spectral function 

In order to study the spectral function in the case when the spin phase fluctuations are taken 
into account, it is convenient to write down the spin functions explicitly in the laboratory system 
of coordinates: 

\a n ) = \a) cos(^ n /2) + 2a\a) sin((^ n /2), (18) 

where tp a is the angle between the directions of the laboratory and local systems of coordinates. 
For simplicity, we assume that the spins lie in the plane of the system and their phase fluctuations 
are small. In this case, one can get the following result for the Green's function: 

G a p{x) = g aP {x)V s P m (x), 

where Q a p(x) is the Green's function of the Nambu spinors and 

V spin {x) = e - ^^)^ ))/ 4 
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is the correlator of the spin phases (for details, see Refs. [[28, 29, 3l|). The phase correlators in the 
last two expressions can be easily obtained from the effective action for the spin phases. Namely, 



D spin (T,r) =exp 



n— — oo 



qdqdtp 1 — cos((jr cos y?) cos(fJ n r) 



(2ir) 



JC - Q 2 



(19) 



where J S pin and IC sp i n are the coefficients in front of the gradient and time derivative terms in 
the effective action Eq. (fT¥|) . One can use these exact expressions to calculate the spin correlation 
functions. However, in order to study the qualitative behavior of the spectral function of the system, 
it is enough to approximate the correlation function by 



t (x) = [9{T* —T) + 9{T - T*) exp(-r/£ spin (T))] , 



-r sp ;„t 



[T I T s pi n ) 



(20) 



where T* is the spin BKT temperature and a sp i ni £ sp m, fspin and T are doping- and temperature- 
dependent parameters. In principle, one also needs to include superconducting fluctuations in 
order to describe the properties of the system at T c < T < Tq. In this case, D spm must be 
multiplied by the corresponding superconducting function D with the space-time dependence 
similar to Eq. (|20|) . We assume that D spm describes the total effect of the spin and superconducting 
fluctuations and use the corresponding parameter notations a, £, r and T instead of ones used in 
Eq. ([2H)l . The values of the parameters can be estimated from experiments (see below). 

In the frequency-momentum representation, the Green's function has the following form: 



G(iw„,k) = -T 



G?q- 



lUJ ri 



T 3 e(q) 



e 2 (q) 



-X> spin (k-q, w„ 



where the Fourier transform of the correlator of the phase fluctuations D spm (k 
be found from Eq. ([20]). In particular, at T < T*: 



D spm {i£l n ,ci)=A 



1 



r 



r q 2 + £-2( T )]a 



r 2 : 



(21) 

•q,w„ -iv rn ) can 
(22) 



where A is a parameter. We shall use this expression to analyze the spectral properties of the 
system. For simplicity, we assume that the inverse time correlation length T is proportional to 
temperature and put T = 0.1T/T, where T ~ T*. 

The spectral properties of the holes can be studied by making the analytical continuation 
iuj n — > lo + irj (jj — > +0) and extracting the imaginary part of the Green's function 



A(u,k) = -(i/7r)ImG(o;,k). 



(23) 



In this case, by using Eqs. pTj) - ((23|) and performing one frequency integration, one can get the 
following expression for the spectral function: 



A(u, k) = A J dq | 



e(g) 



(cc;-S(q)) 2 +r 2 



1 - 



{uj + E{q}) 2 +T 2 
1 



[(k- 



(24) 



where A can be most easily found by using the Green's function zeroth spectral moment sum rule. 
Let us get an approximate expression for the spectral function in the case of = and positive 
frequencies. In this case, the last term in the figure brackets can be neglected, and one gets 



A(uj,k) ~ A{u>) j dip q 



1 



(k 2 - 2ky/2m*(^ + uj)cos(ip q ) + 2m*(fi + uj) + £~ 2 (T) 



(25) 



where 



A(lu) = Am* 



d£[l + sign(£ - n)] 



(w - o 2 + r2 



2m* A 



— + arctan(w/r) 
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As it follows from Eq. ([25]) . the Green's function has a cut-like form, contrary to the Fermi- 
liquid pole-like case. Such a dependence is caused by presence of the spin correlation function 
■pspm (k — q, w n — 0Jm ) i n Eq. (f2Tj) , which smoothes out the (5-function peaks of the spectral function 
that come from the denominator of the fermion Green's function. Therefore, the quasi-particle 
residue Z, which can be defined as the coefficient in front of the Fermi quasi-particle spectral 
function 5-peak, is equal to zero. This means that the system is in a non-Fermi-liquid regime. This 
result is qualitatively similar to the result recently obtained by P.W. Anderson in the case of a 
strongly correlated model [[37] and used to describe the spectral function in BiiSr-2,GaCuiO%+& 
in Ref. [[32]: 

a, - f (u\ sin[(l -p)(tt/2 - tan- 1 ^ - v F k)/T])] 

where f(ui/T) is the Fermi function and p and T are parameters in the effective Green's function: 

G(i,x) = n -. (27) 

|X| — VFt 

In particular, the pre-factor t~ p comes from the contribution of the Gutzwiller projection on the 
single occupied states of the strongly correlated system. It was estimated that p ~ 0.12 and 
r = AT + B{k — hp) 2 in the momentum space, where A and B are constants. 

In fact, in our case we can get an approximate cut-like expression for the spectral function by 
putting (f q = in the expression under the integral in Eq. (|25p and integrating over the angle ip q . 
This can be done since the momentum angle region around <p q = gives the largest contribution 
into the integral. In this case, 



7T 



2 



1 



A(u, k) ~ &irm*A - + arctan(w/r) ^. (28) 



(k - yj2m*{ii + uo)) 2 + £- 2 (T) 



As it follows from this equation, the spectral function has the maximum at u) = £(k) = 
k 2 /(2m*), similar to the free hole case (we put /i = for simplicity), but this function is a smooth 
function, different from the delta-function. Our numerical evaluation of the integral in Eq. (|25p 
show that the approximation Eq. ([28]) is correct only at large a (Fig. [2]). In the case of small a, 
the w[£(k)]-dependence for the spectral function maximum is linear at large £(k). In this case, the 
curve begins at finite value of |k|, which decreases with a increasing. 

From Eq. (|28p . one can find an approximate expression for the spectral function at u> = 0: 

A(. = 0, k) * ^rn*A [{k _ kp? \^ {T)r (29) 
As it follows from this equation, the spectral weight on the Fermi level is defined by £ and a. 

5. Density of states 

The density of states can be obtained from Eq. (HMj) by using the standard expression 



N(u) = J -0^A(uj,k). (30) 

Similar to the previous Section, one can show that the system demonstrates a finite DOS at zero 
frequency and it is defined by the parameters £ and a. 

One can estimate the DOS weight at ui — from Eq. ([29]) : 



d 2 k 1 



N(p = 0) ^ 2 ^/(^ 2[fc2 + e - 2(T)r 

- 2(1-7)^(1-°) [WWeW + lf-*-!]. (31) 
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Figure 2. The spectral function peak frequency tj-dispersion f (k) curve at different values of a. 
This curve approaches the free fermion curve u> = £(k) as a increases. 



From Eqs. (f2"5)) , (f2"5)) and (f3"Tj) one can estimate the values of the phenomenological parameters 
£, a and T by comparing theoretical results with experimental data. In particular, one can get that 
r ~ T and £ is a weakly-dependent function of temperature for a wide temperature range above 
T . 

To conclude this Section, we would like to demonstrate how some of the non-Fermi-liquid 
properties could result from the cut-like structure of the Green's function by using the anomalous 
conductivity as an example. One can roughly estimate the conductivity to be proportional to 
the quasiparticle life-time r(k, w) at |k| = kp. This quantity can be estimated to be inversely 
proportional to the imaginary part of the one-hole self-energy, i.e. r(fcj?,a;) ~ 1/ ImTl(kp, u>). On 
the other hand, since 

JmS(k, u>) 
[ lUJ) ~ (k 2 /2m*-n) 2 + JmE 2 (k, w) ' 

where we have neglected the real part of the self-energy, one immediately gets r(kp, ui) ~ A(kp,ui). 
The frequency dependence of this quantity at different values of a is presented in Fig. [31 where we 
have substracted the frequency-independent part from the spectral function. As it follows from this 
Figure, at low and moderate frequencies the dependence of the conductivity on frequency can be 
approximated by a ~ ui b , where b is a parameter. As it follows from Eq. (f!?5)) , at small frequencies 

A { kp,.)-A(kp,.) ~ _^^J_ tl ,. (32) 

Such a dependence obtained by using a rather rough approximation already indicates that the 
cut-like form of the Green's function can result in a non- Fermi-liquid behavior of the system. 

In order to make a comparison with experiments for conductivity and other quantities in the 
pseudogap phase, one needs to take into account more accurately different properties of the ma- 
terials, like the band structure, the antiferromagnetic spin coupling, which defines J and others. 
Such a comparison is also necessary to the phenomenological parameters for the spin fluctuation 
correlation function Eq. |20)) . 
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Figure 3. The leading frequency-dependent term of the spectral function (in arbitrary units) 
at small u>, k = kp and different values of a. This term is proportional to A(kp, w) — A(kF, 0). 
The frequency uj is given in units of l/(2m*£ 2 ). 



6. Conclusions 

In this paper, we have considered the spectral properties of a phenomenological model of HTSC 
in the underdoped regime by taking into account fluctuations of the phases in the AF spin back- 
ground. Namely, we have considered the temperature evolutions of the spectral function and of the 
density of states. By studying the spectral function, we have shown that its temperature depen- 
dence in the case of HTSCs can be qualitatively described by this model in the case of the proper 
choice of the decoherence time correlation length and other parameters for the spin angle correla- 
tion function. These parameters can be taken from experiments and they are directly connected 
with the microscopic model parameters. Similarly, we have derived and analyzed the expression for 
the density of states. Finally, we have shown that the spin fluctuations can be responsible for the 
anomalous behavior of the conductivity in the underdoped regime. We have compared our result 
for the Green's function with the expression proposed in [El], and have shown that both models 
can describe anomalous properties of underdoped cuprates without using exotic models, like the 
marginal Fermi-liquid etc. The validity of the model studied in this paper to describe all properties 
of cuprates in the pseudogap phase can be tested by taking into account different experimental 
phenomena, what requires a farther investigation. 
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